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MARVIN  MARCUS 

M.  Marcus  has  published  a  total  of  165  research  papers  and  16  books 
and  research  monographs  since  the  beginning  of  his  research  career  in  1954, 
His  initial  work  was  concerned  with  the  behavior  of  solutions  to  systems 
of  ordinary  differential  equations  and  appeared  in  the  Princeton  Annals 
of  Mathematics  Studies,  Contributions  to  Nonlinear  Oscillations  between 
1955  and  1960.  Problems  in  stability  theory  for  such  systems  can  frequently 
be  analyzed  by  studying  the  linear  approximation.  This,  in  turn,  leads 
to  classical  problems  in  eigenvalue  localization.  Most  of  Marcus'  work 
starting  in  the  late  1950's  has  been  in  classical  algebraic  stability  theory 
and  its  supporting  and  derivative  disciplines:  multilinear  algebra;  con¬ 
vexity  theory;  theory  of  inequalities;  group  representation  theory;  nu¬ 
merical  algebra. 

During  the  period  of  the  contract  (October  1,  19 78 -September  30,  1979), 
Marcus'  research  has  been  concerned  with  the  following  general  problem 
which  in  turn  separates  into  two  parts.  The  beginnings  of  this  work  was 
presented  in  an  invited  address  at  Gatlinburg  VII,  Conference  on  Numerical 
Linear  Algebra,  December  11-17,  1977.  Further  developments  were  reported 
on  in  an  invited  address  at  the  New  York  Academy  of  Sciences,  Second  Inter¬ 
national  Conference  on  Combinatorial  Mathematics ,  New  York,  April  4-7,  1978. 
The  most  recent  work  completed  in  the  1978-79  period  was  recently  pre¬ 
sented  at  an  invited  address  at  the  California  Institute  of  Technology. 
Marcus  has  been  invited  to  an  NSF  sponsored  Conference  on  Linear  Algebra 
at  Auburn  to  present  a  survey  lecture  on  this  research  and  later  develop¬ 
ments  (March  1980). 
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GENERAL  PROBLEM.  Let  V  be  an  n-dlmensional  unitary  space  with  inner 
produc t  ( •,  •)  and  let  T:V  -*  V  be  a  linear  operator.  Given  a  function 
f:C  “•  C  ,  investigate  W^(T)  ,  the  set  of  points  in  the  complex  plane,  C  , 
of  the  form 

f  ( (Tx^,x^) p  • .  •  y  (Tx^,x^) ^  • 

The  vectors  x^,  . . . ,  x^  usually  run  over  prescribed  sets  of  k  orthonormal 
vectors  in  V  ;  however,  sometimes  the  hypotheses  stipulate  the  value  of  the 
inner  product  (x^,Xj)  to  be  something  other  than  6^  ,  i, j  *  1,  . .  .,  k  . 

If  k  ■  1  and  f(z)  ■  z  then  W(T)  ,  the  classical  Toeplitz-Hausdorf f 
numerical  range  of  T  ,  is  a  specific  instance  of  W^(T)  .  Thus  W^(T) 
will  be  called  the  f -numerical  range  of  T  . 

This  General  Problem  separates  into  the  following  specific  types  of 


investigations : 

I.  The  relationship  between  the 


of  W^(T)  and  the  eigenva lues , 


singular  values ,  elementary  divisors ,  etc .  of  T  . 
II.  Numerical  computations  of  W-(T)  and  related  sets . 


Results  related  to  I. 

In  a  recent  query  A.  Abian  posed  the  following  interesting  question. 
Let  V  be  an  n-dimensional  inner  product  space,  and  let  A  ,  B  ,  P  ,  Q  be 
linear  on  V  .  What  relations  must  exist  among  these  operators  so  that 
the  inequality 

(Av,u) (Bu,v)  £  (Pu,u)(Qv,v)  (1) 


holds  for  all  u  and  V  ? 


Marcus  solved  this  problem  completely  in  [7]. 
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THEOREM*  Let  V  be  a  unitary  space,  dim  V  ^  3  ,  and  assume  that  A,  B,  P,Q 
are  nons ingulag .  Then  (1)  holds  for  all  u  and  v  if  and  only  if 

(i)  P  -  OH  ,  Q  *  0K  ,  aP*  e  ■  ±1  ,  H  and  K  are  definite  hennitian  and 
(ii)  A*  -  Xb ,  X  real,  so  that  (1)  reads 

X|(Bu,v)|2  3  e(Hu,u)(Kv,v)  ,  (2) 

and 

(iii)  if_  e  =  1  ,  X  >  0  ,  then  H  and  K  have  the  same  sign  (i.*e., 
both  positive  definite  or  both  negative  definite)  and 

X  (p”1AQ'1B)  5  1  ;  (3) 

max 

or 


(iv) 

If 

e  = 

i ,  x  < 

o, 

then 

H 

and 

K 

have 

the  same  sign ; 

or 

(v) 

If 

e  = 

-1  ,  X  > 

o. 

then 

H 

and 

K 

have 

opposite  signs 

and 

X  (P-1AQ-1 

max 

B)  S 

i  ; 

w 

(vl) 

if 

e  30 

f 

Vi 

A 

0 , 

then 

H 

and 

K 

have  opposite  signs 

This  result  was  proved  using  properties  of  linear  maps  on  the  tensor 
space  over  V.  The  same  techniques  led  to  results  of  the  following  type: 

Let  A p  be  linear  maps  on  V,  dim  V  ■  n  ,  and  define  the 

set 

k 

W(A1,...,Ak)  ][(AiVx1)|V...,xk  o.n.j  . 

1-1 

THEOREM.  IT  A^,  . . . ,  A^  are  linear  maps  on  a  unitary  space  V  ,  k  5  n  - 
dim  V  ,  then 

W<A1,...,Ak)  -  (0) 
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iff  some  i£  0  [ 12} . 

Using  properties  of  the  Grassmann  algebra  the  following  results  were 
obtained  in  [10] . 

Let  A  be  an  n  X  n  normal  matrix  over  C  ,  and  9^  n  be  the  set 
of  strictly  Increasing  Integer  sequences  of  length  m  chosen  from 
1,  . n  .  For  a,  P  ^  Sn  n  denote  by  A[a|P]  the  submatrix  obtained 
from  A  by  using  rows  numbered  Ot  and  columns  numbered  p  .  For 
k  €  {0,1,... ,m)  write  |afip|  M  k  if  there  exists  a  rearrangement  of 
1,  ...,  m  say  iy  ...,  ,  i^,  ...,  i^  such  that  0(1^)  *  0(ij)  , 

i  -  1,  •*.,  k,  and  (0(1^),  .  .  .,<*(im)  }  H  {PUk+1),  •  •  •  )  -  0  .  A 

new  bound  for  |det  A[a|P]|  is  obtained  in  terms  of  the  eigenvalues  of 
A  when  2m  =  n  ,  and  |  or  0  0  j  *  0  . 

THEOREM.  Let  m  >  2  ,  n  *  2m  ,  and  let  A  be  an  n  X  n  normal  matrix 

over  C  with  eigenvalues  X^ ,  . . , ,  X^  .  Denote  by  X^  the  product 

X  •  •  •  X  ,  N  ,  for  0)  €  Q  .  Then 
ij(l)  a(m)  '  -  m,n  - 


|det  A[U)|ou# ] | 


i  fE JW 

uu  €Q 


T~T~j,  .  V  I  |\  +  (-1)"V  «| 

2  (m+1)  w  €  1  uu  »  1 


if  m 


if  tn  >  2 


The  set  Q  c  ^  is  any  subset  with  j Q !  3  ~(“)  • 

For  example,  if  m  *  2  ,  n  *  4  take  Q  -  { (12),  (13),  (14) }  ,  Q*  = 

{ (34), (24) , (23) ]  and  the  result  states  that  jdet  A[ 12 | 34 ] |  ^ 


JiiWVJ  + 


|  ^1^2  +  ^“2^4  I  +l  ^1^4  *  ^2^3  I  ^ 


In  [4]  the  following  results  are  obtained.  Let  E  (|X|)  = 


m 


E^( | Xi | , . . ., | Xn | )  denote  the  mth  elementary  symmetric  polynomial  in 


i.e., 
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n 

Em(lXl)  9  L  IT  I  Vi)  *  * 

0)€Q  1*1 

m,  n 

THEOREM,  Let  A  be  a  normal  matrix  with  eigenvalues  X^ ,  . . . ,  X^  ,  n  £  4 . 
Let  B  be  an  m-square  submatrix  of  A  having  precisely  k  main  diagonal 
entries  lying  on  the  main  diagonal  of  A.  lf_  k^m-2  then 

|det  B|  S 


Some  corollaries  that  are  consequences  of  this  result  are  directly 
applicable  to  eigenvalue  localization  problems. 

Recall  that  the  spread  of  A  is  the  number 

s  (A)  *  max  I  X  -  X.  I  . 

i,j  1  ] 

Corollary  1.  Lf  2  <  m  <  n  and  the  m- square  submatrix  B  of  A  has 
no  main  diagonal  elements  lying  on  the  main  diagonal  of  A  then 

8(A)  £73^2(0*1)^)  Y/m  |det  B|1/m  ,  (5) 

and  if  m  =  2 

s (A)  S  2v^(n(n-l))‘1/2  |det  B|1/2  . 

In  the  second  corollary  A  is  an  arbitrary  n-square  matrix,  (n^2)  . 

Let  d  be  the  largest  m-square  subdeterminant  of  A  (in  absolute  value) 
m 

and  let  ||a|1  be  the  Hilbert  norm  of  A,  l.e.,  the  largest  singular  value 


of  A  . 
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Corollary  2.  I_f  2  <  m  £  n  then 

and 

ltA||  &  2(n(n-l))”l/2  d\'2  . 

Results  related  to  II 

In  [5]  the  following  results  were  obtained: 

THEOREM.  Let  A  be  a  0, 1 -matrix  with  at  most  one  1  in  each  row  and 
column.  Then  the  numerical  range  of  A  jLs_  the  convex  hull  of  a  polygon 
and  a  circular  disk,  both  centered  at  the  origin  and  contained  in  the  unit 
disk. 

The  proof  uses  a  permutation  similarity  to  reduce  A  to  a  direct 
sum  of  matrices  whose  numerical  ranges  can  be  determined.  A  computer  pro¬ 
gram  was  developed  which  plots  the  boundary  of  the  numerical  range  of  an 
arbitrary  complex  matrix.  For  example,  for  n  -  10  consider  the  matrix 


A  - 


000000000 


100000000 


000100000 


001000000 


000000000 


010000000 


000000000 


000000100 


000010000 


0000000010 
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To  exhibit  the  sensitivity  of  W(A)  to  changes  in  A  let  B  be  the 
matrix  obtained  from  A  by  replacing  the  0  in  the  (5,6)  position  by  a  1. 
The  following  are  plots  of  W(A)  and  W(B)  produced  by  a  program  using 
standard  FORTRAN  subroutines  and  run  on  the  AS/6  at  the  UCSB  Computer  Center. 


An  unresolved  conjecture  of  Marcus  [10]  states  that  for  a  normal  matrix 
with  given  eigenvalues  the  closer  a  subdeterminant  is  to  being  principal, 
the  larger  its  absolute  value  can  be.  More  precisely  let  A  be  an  nxn 
matrix  and  consider  the  m  X  m  subdeterminants  of  U*AU  lying  in  rows 
1,  . ..,  m,  columns  1,  . ..,  k  ,  nrt-1,  . 2m-k  : 

jdet  U*AU[1, . • . ,m | 1, . . • ,k,m+l, . . • ,2m-k] |  ; 

U  runs  over  all  unitary  matrices.  The  columns  labeled  k  =  0  ,  k  *  1  , 
k  *  2  ,  k  =  3  and  k  =  4  in  the  following  Table  contain  the  maximum 
values  obtained  for  these  numbers.  One  hundred  random  unitary  matrices 
were  generated  for  each  run:  The  eigenvalues  of  A  (the  only  pertinent 
information)  are  listed  in  the  left  hand  column.  The  100  ‘random1  uni¬ 
tary  matrices  were  generated  as  follows: 
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1.  n(n+l)/2  random  complex  numbers  were  generated  to  obtain  an 
n  X  n  skew-herroltlan  matrix  S  . 

2.  Then  standard  FORTRAN  subroutines  were  used  to  compute  the  Cayley 
transform 

U  =  (in- S)(In  +  S)_1  , 

a  unitary  matrix. 
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PUBLICATIONS 

completed  under  AFOSR  Contract  No.  AF  F4962078-C-0030 
October  1,  1978  -  September  30,  1979 
MORRIS  NEWMAN 

1.  Combinatorial  matrices  with  small  determinants,  Canad.  J.  Math.  30 
(1978),  756-762. 

2.  Normalizers  of  modular  groups.  Math.  Ann.  238(1978),  123-129. 

3.  Equivalence  without  determinantal  divisors,  Linear  and  Multilinear 
Algebra  7(1979),  107-109. 

4.  The  use  of  integral  operators  in  number  theory  (with  C.  Ryavec  and 
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11.  On  a  problem  suggested  by  Olga  Taussky-Todd,  to  appear  in  Illinois 
J.  Math. 

12.  Positive  definite  matrices  and  Catalan  Numbers  (with  F.T.  Leighton), 
to  appear  in  Proc.  Amer.  Math.  Soc . 

13.  Matrices  of  finite  period  and  some  special  linear  equations,  to  appear 
in  Linear  and  Multilinear  Algebra. 

14.  Gersgorin  revisited,  to  appear  in  Letters  in  Linear  Algebra. 

15.  A  surprising  determinantal  inequality  for  real  matrices  (with  C.R. 
Johnson),  to  appear  in  Math.  Ann. 
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17.  Determinants  of  circulants  of  prime  power  order,  to  appear  in  Linear 
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MORRIS  NEWMAN 

The  principal  object  of  Newman^  research  is  to  apply  number- theoretic 
ideas  to  problems  of  numerical  analysis  and  computation.  In  particular, 
applications  of  modular  arithmetic  are  being  made  to  such  topics  as  (1) 
the  exact  solution  of  an  integral  system  of  linear  equations,  (2)  the  de¬ 
termination  of  the  exact  inverse  of  an  integral  matrix,  (3)  the  determi- 
nation  of  the  rank  and  a  basis  for  the  null  space  of  an  integral  matrix, 

(4)  the  determination  of  the  eigenvalues  of  a  rational  triple  diagonal 
matrix  to  any  desired  degree  of  accuracy,  and  (5)  the  computation  of  the 
permanent  of  a  matrix.  Programs  to  perform  (1)  and  (2)  which  take  full 
advantage  of  parallel  computation  have  been  prepared  for  ILLIAC  IV.  For 
example,  the  exact  solution  of  any  integral  linear  system  in  at  most  40 
unknowns,  together  with  the  determinant,  can  be  found  in  less  than  10 
seconds.  ILLIAC  IV  programs  for  (3),  (4),  and  (5)  are  in  process  of  pre¬ 
paration  . 

Past  Accomplishments  Under  The  Existing  Contract 

The  primary  objective  of  the  research  undertaken  in  the  last  period 
was  to  produce  an  ILLIAC  IV  program  to  find  the  exact  inverse  of  an 

n  X  n  nonsingular  matrix  A  with  integer  entries.  This  has  been  achieved. 

2 

The  program  requires  just  n  + 2n  cells  in  each  processing  element  of 
ILLIAC  IV  for  storage,  and  takes  full  advantage  of  the  parallel  processing 
capabilities  of  the  machine.  For  example,  the  exact  inverse  of  any 
40  X  40  matrix  of  this  type  can  be  found  in  approximately  30  seconds. 

The  program  actually  produces  the  determinant  and  adjugate  matrix,  and 
may  be  combined  with  another  existing  program  to  compute  A~^B,  where  B 
is  any  n  X  k  integral  matrix.  The  time  to  perform  a  similar  calculation 
on  a  serial  machine  could  well  be  as  high  as  3  hours. 


Con  feta  with  Air  Force  Laboratories 

Newman  has  had  invaluable  help  from  the  staff  of  the  Institute  for 
Advanced  Computation  at  Sunnyvale,  and  has  also  discussed  the  project 
with  Paul  Nikolai,  at  Wright-Patterson  Air  Force  Base. 
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Diagonal  Elements  of  Complex  Symmetric  Matrices 

It  was  proposed  to  study  the  relations  between  the  diagonal  elements 

and  singulars  of  complex  symmetric  matrices,  (Recall  that  the  singular 
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values  of  a  matrix  A  are  the  eigenvalues  of  (AA*)  ,  A*  the  conju¬ 

gate  transpose  of  A  .  )  This  problem  Arose  from  elementary  particle 
physics ,  where  the  matrix  A  in  Fermi  statistics  is  complex  synraetric 
and  unitary.  A  complete  solution  was  obtained,  the  form  of  the  solution 
being  unusually  striking  and  unexpected.  Specifically,  let  the  diagonal 
elements  be  d^  ,  d2>  .  d^  ,  numbered  such  that  |d^|  ^  £  |d^|  , 

and  let  the  singular  values  be  8j  '  ^  s^  .  (These  are  nonnegative 

numbers.)  Then  a  complex  symmetric  matrix  exists  with  diagonal  elements 
d^,  d^  and  singular  values  s^,  . ..,  s^  ,  if  and  only  if 
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and  one  further  singularity. 
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(The  last  condition  has  three  minus  terms  on  the  left  and  two  on  the  right.) 

The  Important  part  about  this  result  is,  not  so  much  its  unexpected 
form,  but  the  fact  that  connections  with  the  structure  theory  of  simple 
Lie  algebras  are  suggested.  It  is  plain  that  underlying  these  Inequalities 
is  a  root  system  for  some  simple  Lie  algebra,  and  work  is  now  under  way 
to  identify  this  connection. 

Matrix  Valued  Inequalities  and  Equalities;  an  Exponential  Formula 

During  the  term  of  the  grant,  the  proposer  was  led  to  the  following 
conjecture. 

Let  /  and  B  be  matrices,  real  or  complex,  of  the  same  size. 

Then  one  would  like  to  have 

A  B  A+B 
e  e  *  e 

But  of  course  this  rarely  happens,  unless  the  matrices  commute.  However, 
on  the  basis  of  his  experience  with  matrix  valued  inequalities,  (see  below), 
the  proposer  conjectured  that 

A  B  UAU*  +  VBV* 
e  e  3  e  , 

for  suitable  unitary  matrices  dependent  on  A  and  B  .  Considerable 
work  has  been  carried  out  on  this  conjecture,  with  these  results: 

(a)  the  conjecture  has  been  verified  in  many  cases, 

(b)  it  has  been  formally  verified,  meaning  that  infinite  series  for 
U  and  V  have  been  found  that  work  when  used  in  conjunction  with  the  in- 
finite  series  for  e  • 

All  that  remains  is  to  prove  that  the  infinite  series  giving  U  and 
V  are  convergent.  They  seem  to  converge  in  examples,  but  strenuous  ef¬ 
forts  to  establish  this  in  general  have  so  far  not  succeeded. 
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This  problem  is  basically  Lie -theoretic,  the  central  tool  being  the  Campbell 
Baker  Hausdorff  formula  from  Lie  theory. 

The  Matrix  Valued  Triangle  Inequality 

1/2 

If  A  Is  a  matrix,  define  its  matrix  absolute  value  as  j A J  *  (AA*) 

Let  B  be  another  matrix.  In  analogy  with  scalars,  one  %*>uld  like  to 
have  the  triangle  inequality 

|A+B|  S  |A|  +  |  B|  , 

meaning  that  the  right  side  minus  the  left  side  is  positive  semidef inite . 
This  inequality  is  usually  false.  Hie  proposer  has  shown  that  it  becomes 
true  if  it  is  changed  to 

|A  +  B|  5  u  |  A  |  u*  +  v|b|v*  (1) 

where  U  ,  V  are  unitary  matrices  dependent  on  A  and  B  .  During  the 
course  of  the  grant,  the  proposer  discovered  that  this  inequality  shares 
all  the  properties  of  the  scalar  triangle  inequality: 

Equality  holds  in  (1)  if  and  only  if  A  and  B 
have  polar  factorizations  with  a  conmon  unitary  factor. 
Further  work  is  continuing  on  matrix  valued  inequalities. 

Other  Work 

A  variety  of  related  problems  involving  eigenvalues,  singular  values, 
matrix  inequalities,  have  been  examined  and  continue  to  be  examined.  There 
are  ongoing  efforts  to  tie  up  these  investigations  with  Lie  theory,  with 
some  success  to  date,  but  probably  with  more  success  in  the  future. 
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1.  On  a  conjecture  of  R.  F.  Scott  (1881),  Linear  Algebra  and  Appl.  28 
(1979),  141-153. 

2.  An  asymptotic  solution  of  the  multidimensional  dimer  problem,  Linear 
and  Multilinear  Algebra,  to  appear. 

3.  Bounds  for  permanents  and  determinants,  Linear  and  Multilinear  Algebra, 
to  appear. 

4.  Rearrangement  inequalities,  to  appear. 


HENRYK  MINC 


A  substantial  part  of  Mine's  research  concerns  permanents  and  their 
applications.  In  paper  [1]  an  old,  outstanding  conjecture  of  R.F.  Scott 
is  proved.  This  is  a  formula  for  the  permanent  of  a  special  type  of 
Cauchy  matrices.  It  is  proved  by  evaluating  the  eigenvalues  of  certain 
circulants  and  their  determinants.  These  results  are  of  independent  in¬ 
terest.  In  [2]  a  new  lower  bound  for  the  dimer  problem  is  obtained.  Al¬ 
though  the  bound,  together  with  an  upper  bound  previously  obtained  by  Mine, 
leads  to  an  asymptotic  solution  of  the  dimer  problem,  the  all-important 
three  dimensional  case  remains  unsolved.  In  order  to  improve  the  known 
bounds  for  the  three-dimensional  problem,  a  large  number  of  permanents  of 
(0, l)-circulants  were  computed  and  tabulated,  but  conclusive  results  have 
not  been  obtained  yet. 


Schinzel  obtained  the  following  remarkable  inequality  for  the  deter¬ 


minant  of  a  real  matrix 


A  = 


(a 


|  det (A) 


-  TM  X 


q  =  l 
a  >0 

pq 


n 


q~i 

a  *0 

pq 


Schinzel's  formula  was  substantially  improved  by  Johnson  and  Newman.  In 
[3]  their  formula  was  used  to  obtain  bounds  for  the  absolute  value  of  the 
determinants  of  complex  matrices.  Johnson  and  Newman  also  applied  their 
method  to  permanents  of  real  matrices  and  showed  that  if  A  is  an  n  X  n 
matrix  such  that  the  sum  of  the  positive  entries  and  the  absolute  value 
of  the  negative  entries  m  each  row  do  not  exceed  1  ,  then 


| per(A) |  S  2 


I 
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In  [3]  a  rearrangement  inequality  is  obtained,  and  used  to  strengthen 
the  bound  in  (*) .  This  rearrangement  inequality  is  generalized  in  [4] 


where  eight  other  new  rearrangement  inequalities  are  obtained.  For  ex¬ 


ample,  it  is  shown  that  if  (a  ,  ...,a  , )  are  nonnegative  m  -tuples, 

Ij  ®j  > J  J 


* ,  n  ,  m.  £  •  •  ♦  £  m  &  1  ,  then 
1  n  9 


This  result  and  the  other  inequalities  in  [4]  generalize  and  extend  re¬ 
arrangement  inequalities  obtained  by  Mine  in  1971. 
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Interscience,  1978,  Linear  and  Multilinear  Algebra  6(1978),  331-338. 

2.  Book  Review:  On  Systems  Analysis:  An  Essay  Concerning  the  Limitations 
of  Some  Mathematical  Methods  in  the  Social,  Political  and  Biological 
Sciences ,  David  Berlinski,  MIT  Press,  1976,  Linear  and  Multilinear 
Algebra  6(1978)  ,  333-335. 

3.  Book  Review:  Applied  Linear  Algebra,  Second  Edition,  Ben  Noble  and 
James  Daniel,  Prentice  Hall,  Linear  and  Multilinear  Algebra  6(1978), 
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4.  Inversion  of  order  in  non-absolutely  convergent  multiple  integrals, 
to  appear. 

5.  Solution  of  the  Kac-Siegert  equations,  to  appear. 

6.  Application  of  the  Kac-Siegert  equations  to  distributions  of  spectral 
estimates,  to  appear. 

7.  Sampling  versus  quantizing  in  stochastic  processes,  to  appear. 

8.  Schwarz  distribution  theory  for  linear  difference-differential  equations, 
to  appear. 
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The  scalar  discrete  optimal  gain  solution  discussed  in  the  1977-78 
annual  report  was  extended  to  2x2  matrices  with  diagonal  objective 
functions.  Linear  control  of  non-linear  systems  was  accomplished  for  the 
case  of  2x2  quadratic  systems  of  Hille  (prey-predator)  type  -  this  has 
ecological  applica tions .  Non-linear  control  of  3x3  linear  systems  was 
accomplished  for  the  case  of  conflict  between  unsymmetr ical  adversaries 
(one  aggressive,  one  defensive)  and  simulated  for  symmetrical  adversaries 
(both  aggressive).  The  linear  control  of  space-craft  attitude  (non-linear 
equations)  was  illuminated  by  numerous  analog  computer  studies  but  a  satis¬ 
factory  theory  remains  elusive. 

The  control  of  saturating  lasers  was  studied  by  an  integral  equation 
method  of  T.  Triffet  (of  Arizona, Tucson)  and  H.S.  Green  (of  Adelaide, 
Australia).  The  results  can  be  applied  to  genetic  problems  in  which  the 
reproduction  ratios  are  variable,  decreasing  (but  not  to  zero)  as  population 
densities  increase.  The  details  involve  the  spectra  of  rather  complicated 
differential  operators-in  general,  it  can  be  said  that  the  saturation  is 
not  exponential,  but  can  be  approximated  by  a  combination  of  several  ex¬ 
ponentials;  the  decay  coefficients  are  expressible  in  terms  of  physically 
(or  genetically)  significant  quantities. 

Some  attention  was  given  to  obtaining  the  characteristic  function  of 
the  lognormal  distribution,  and  to  finding  the  distribution  of  the  sum  of 
independent,  but  not  identical,  lognormal  random  variables.  This  work 
has  numerous  applications  in  radar  detection  and  component  reliability. 

The  results  are  expressible  in  terms  of  series  of  Hermite  functions  with 
gamma  (or  incomplete  gamma)  function  coefficients. 


